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. Abstract 

, New solvable vertex models can be easily obtained by staggering the spectral parameter in already known 

ones. This simple construction reveals some surprises: for appropriate values of the staggering, highly non-trivial 
continuum limits can be obtained. The simplest case of staggering with period two (the Z2 case) for the six- 
vertex model was shown to be related, in one regime of the spectral parameter, to the critical antiferromagnetic 
Potts model on the square lattice, and has a non-compact continuum limit. Here, we study the other regime: in 
^ ' the very anisotropic limit, it can be viewed as a zig-zag spin chain with spin anisotropy, or as an anyonic chain 

m : with a generic (non-integer) number of species. From the Bethe-Ansatz solution, we obtain the central charge 

c = 2, the conformal spectrum, and the continuum partition function, corresponding to one free boson and two 
Majorana fermions. Finally, we obtain a massive integrable deformation of the model on the lattice. Interestingly, 

its scattering theory is a massive version of the one for the flow between minimal models. The corresponding field 

(2) 

theory is argued to be a complex version of the C2 Toda theory. 
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• 1 Introduction 



It is a simple consequence of the quantum inverse scattering [Tj formalism (going back to Baxter's 'Z invariance' [5]) 
^ . , that new integrable vertex models can be obtained from basic ones by allowing for some staggering of the spectral 
' parameters. If the basic i?-matrix is associated with a single crossing , one can in this way build 'block' _R-matrices, 
using crossings, with n = 2,3,..., and staggering the spectral parameters (see Fig. [T|). 

While constructing the model and writing down the Bethc equations is straightforward, the physics of the models 
thus obtained presents interesting subtleties. A striking example was discussed in detail in Rcfs. |31I3], in relation 
with the antiferromagnetic Potts model [5] . It was found then that the case n = 2 for the six- vertex model has a non- 
compact continuum limit [3] in a certain regime of the spectral parameter (see below for a more accurate definition), 
and may be related to the complex sine-Gordon (SG) model. While major difficulties remain in this case, the other 
regime of the spectral parameter turns out to be also of interest, and somewhat more tractable. Its study is the main 
goal of this paper. 

Before launching into details, a general discussion about vertex models, spin chains and field theories seems useful. 
Indeed, the correspondence between integrable spin chains with SU(2) or SU(2)q symmetry and quantum field theories 
has been investigated in great details already. In the antiferromagnetic regime, it is well known that a chain of spin s 
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Figure 1: The 7?-niatrices for a basic vertex model (left), and for a staggered vertex model with n = 2 (right). 



corresponds to a lcvel-2s Weiss- Zumino-Wittcn (WZW) model in the SU(2) case, and to a deformation of this theory 
in the Cartan direction for SU(2)(j [B]. Cursory examination of the literature would suggest nothing much remains to 
be done in this area. 

From the i?-matrix point of view, models with higher values of the spin s are obtained by projecting tensor products 
of (2s) fundamental representations of U(j(sl2) onto higher spin. This construction has a close parallel in conformal 
field theory (CFT), where higher level (2s) representations of current algebra are obtained by combining (2s) level 
one representations. There are several reasons why it would be interesting to build integrable models which are not 
projected onto irreducible Ug(sl2) components. In the case s = 1 for instance, this would correspond to models of pairs 
of spin-i variables. This can be reinterpreted more physically in terms of ladders, or, in terms of electron physics, 
pairs of wires or channels. The latter case is of the highest importance. For instance, the two-channel Kondo model [7] 
or the two-channel Interacting Resonant Level models [H [9] are usually solved by going to an even-odd basis, which 
effectively amounts to solving the problem in the level k = 2 sector. Many physical questions are however related to 
the mixture of the even and odd degrees of freedom - e.g. because it corresponds to transport of electrons between 
wires. The search for integrable cases where this mixture could be studied is a priority. 

To have a better idea of what to expect, it is useful to turn to the CFT point of view. Imagine starting with two 



SU(2) Kac-Moody algebras at level 1, represented by the currents j[ 



1, 2 and ji — 1,2, 3. The sums 



(1.1) 



are well known to provide then a Kac-Moody algebra at level 2. Of course, each level 1 corresponds to central charge 
c = 1, while level 2 has central charge c = 3/2. The point is that in taking two copies of level 1, an Ising model CFT 
factors out, according to the well known-decomposition 



SU(2)i X SU(2)i = SU2(2) + Ising. 

This can be illustrated quickly using bosonisation. Introduce two chiral bosons 0i , 02 with propagators 
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The two level- 1 current algebras are obtained through 
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It is convenient to introduce now symmetric and antisymmetric combinations of the bosons: 

^(01+02), = (01 - 02) ■ 



(1.2) 



(1.3) 



(1.4) 



(1.5) 
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So we have 

J oc cos V4^(/) exp(±i%/4^$) , oc 9$ . (1.6) 

The field cos(-\/47r'i>) is a Major ana fermion |TU]. The field i sin(-\/47r(^) is another one, which is orthogonal to the 
currents J'^, and is discarded in the construction of SU(2)2. Corresponding to this splitting, the sum H = Hi + H2 oi 
the two one-boson Hamihonians decomposes as H = -ffsu(2)2 + -Rising where i?su(2)2 ^ J2fj, ■ J^J^ ■, and are the 
currents at level two. 

Quantum deformations of SU(2)2 arc obtained by adding to the Hamiltonian i?su(2)2 ^ Cartan deformation pro- 
portional to J^J^. We can as well deform the full Hamiltonian, obtaining in this way a theory made of two Majorana 
fermions and one boson with anisotropy-dependent radius. This should be the continuum limit of the models we are 
after. We will show in the following that these models are obtained by the general staggering construction, with n = 2 
and an appropriate choice of the spectral parameters. 

Interestingly, it can be shown [11] that the staggered models correspond algebraically to solutions of the Yang- 
Baxter equations based on 'bigger' irreducible representations of the quantum affine algebra Uq(sl2). This occurs 
ultimately because finite-dimensional irreducible representations of quantum affine algebras are isomorphic to products 
of evaluation representations, which arc themselves 'decorations' (with the spectral parameter) of the usual spin-s 
representations of Ug(sl2) [T2j . 

Another important application of the staggered model is that it can be used to produce a lattice discretisation 
of a massive QFT. This is done, following Ref. [13], by introducing into the staggered model an additional (purely 
imaginary) staggering of the spectral parameters. Using this approach, we obtain a scattering theory involving two 
types of massive particles, where the scattering between particles of the same type (resp. different types) is given by 
the Sine-Gordon S'-matrix (resp. the Sine-Gordon S'-matrix with an imaginary shift in the rapidity). It turns out 
that such a scattering theory (but with massless particles) arose before [14] in the context of minimal models of CFT 
perturbed by the $13 operator. Consider the action: 

A = Anun + xJd^X'i>i3{x), (1.7) 

where Amin is the action of a minimal model of CFT. The problem of finding the Thermodynamic Bcthe Ansatz (TEA) 
equations for the renormalisation-group fiow of the theory (|1.7p was first studied by Zamolodchikov in Ref. [TS] [16] . 
The results depend crucially on the sign of the coupling A. For A < 0, the model (|1.7p becomes massive. It was 
shown in Ref. [15j that the corresponding S'-matrix is the simple RSOS S'-matrix and that the TBA diagram is of 
the An type, with a massive particle at one end of the diagram. For A > 0, the model p.7p describes the massless 
fiow between two consecutive minimal models. In Ref. [TB], a TBA diagram was proposed, without resorting to an 
S'-matrix: this diagram is also of the An type, but with mass terms e^* at the two ends of the diagram. It was found 
later, in Ref. [ij, that the corresponding scattering theory consists in massless left/right particles, interacting through 
the SG and shifted-SG S'-matrices. 

In summary, we construct here a non-critical lattice model whose continuum excitations are described by a massive 
version of the S-matrix for the flow between minimal models of CFT. We also propose an effective QFT for this model, 
with one boson and two Majorana fermions which interact with each other. 

The paper is organised as follows. In Section [2] we expose in more detail the construction of the model and its 
various lattice formulations, as well as the relation with spin chains and anyonic chains. In Section [31 we present the 
Bethe- Ansatz solution, and obtain the critical exponents, through the study of low-energy excitations. In Section [U 
we discuss the associated CFT, and exhibit the full operator content through the study of torus partition functions. 
Finally, in Section [5] we study the integrable massive deformation, and derive its scattering theory, TBA equations 
and ground-state energy scaling function. This allows us to propose an interacting effective QFT. Some important 
but quite long calculations are done in the Appendices. 
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Figure 2: Equivalence between the cluster and loop models. The full dots represent Potts spins, and the empty dots 
are the sites of the dual lattice. The thin lines arc the edges of the lattice £. Each cluster configuration corresponds 
to precisely one loop configuration. 



2 Solvable Potts models and loop Hamiltonians 

In this section, we first recall the definition of the Potts model [17] on the square lattice and its equivalence to a loop 
model based on the Temperley-Lieb (TL) algebra HH]. Then we introduce a solvable case, involving a staggering of 
the spectral parameters, and we obtain the expression (|2.18|) for the associated one-dimensional Hamiltonian, in terms 
of the TL generators. This solvable model has two critical regimes: regime I corresponds to the antiferromagnetic 
critical transition [31I11[S]; regime II contains a mixture of ferromagnetic and antiferromagnetic couplings, and is the 
subject of the present paper. Finally, we explain the relation to Majumdar-Ghosh spin chains ^IM [20] and anyonic 
chains [HJUg. 

2.1 The Potts model and the Temperley-Lieb algebra 

The Q-statc Potts model [TJ is a model of classical spins with nearest-neighbour interactions. Each spin Sj can take 
Q values, and sits on a vertex of the square lattice. The Boltzmann weights and the partition function are given by: 

W[{S,}] = l[cxp[J,^ S{S,,Sj)], 

(2.1) 

Zpotts = il- 
ls,} 

where the product runs over all the bonds of the lattice, and the Jij are the coupling constants of the model. This 
model can be reformulated as a cluster model, called the Fortuin-Kasteleyn model [53] , in the following way. We may 
write the Boltzmann weights as: 

W[{S,}] = n [1 + ^^ASr,S,)] , V,, EE e-'-- - 1 . (2.2) 

iv) 

The above expression can be expanded, and each term in the expansion is represented by a subgraph of the square 
lattice, consisting of the bonds where the 6 term has been chosen. When we sum over the spin configurations, the 
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Figure 3: The Temper ley-Lieb algebraic relations (|2.6 





subgraph gets a weight Q for each of its connected components. So the partition function can be written as: 



^FK 



Q"^"' n ^^^ ^ 



G 



(2.3) 



where the sum is over all possible subgraphs (or cluster configurations) G, and C{G) is the number of connected 
components of G. The Fortuin-Kasteleyn model is, in turn, transformed to a dense loop model. Consider the lattice 
C obtained by the union of the original square lattice and its dual. Each face of C contains exactly one Potts bond. 
By decorating these faces as shown in Fig. [2l we obtain configurations of closed loops on the lattice. The partition 
function is then: 



where L{G) is the number of closed loops. 

In the present study, we will restrict ourselves to the case when the couplings are Jy 



Ji on horizontal bonds 



and Jij = J2 on vertical bonds. Then, for the loops, a plaquette can be of two types, according to the direction of the 
Potts bond it contains. The local Boltzmann weights for the loop plaquettes can be written: 



(2.5) 




where x{u) = vi/^/Q for plaquettes of type 1 and x{u) - 
closed loops are encoded in the Temperley-Lieb algebra 

'Q e,j 



Q/v2 for plaquettes of type 2. The non-local weights for 
with generators {ej)j=i...2N' 



(2.6) 



ejef=ej,ej if|j-j'|>l 



These relations are interpreted graphically as follows: j is the site index in the horizontal direction, and Cj is the 
represented by the second term of Eq. (|2.5p at position j (see Fig. [3]). The loop weight ^/Q is parametcrised as: 



Q = 2 cos 7 , 



0<7<2 



(2.7) 
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Figure 4: The Yang-Baxter equations (|2.8p . 



2.2 Solvable inhomogeneous Potts model 

To construct a solvable model, the first step is to obtain an 7?-niatrix which satisfies the Yang-Baxter equations: 

Rj,j+iiu)Rj+i,j+2iu ~ v)Rj^j+i{v) = Rj+i^j+2iv)Rj,j+iiu - v)Rj+i,j+2{u) . (2.8) 
The Temperley-Lieb algebra (|2.6p provides a solution to these equations: 

Rj_j+i{u) = sin(7 — u) 1 + sinu Cj . (2-9) 
Comparing with (|2.5[) , the local weight x is given by: 

<u)= . (2.10) 
sm(7 — u) 

Using this _R-matrix, we can construct a solvable model on the square lattice, by choosing the spectral parameters 
along the lines of the lattice. Suppose, from now on, that the square lattice where the loops live consists of horizontal 
and vertical lines. To respect the alternation Ji, J2 of the Potts coupling constants, we have to use the spectral 
parameters {u,u + a^u,u + . . .) and (0, a, 0, a, . . . ), and also ensure that R{u -I- a) cx R{u — a). This holds only for 
a = or a = 7r/2. 

The case a = is the homogeneous TL loop model, and corresponds to the well-studied self-dual Potts model. 
It has two critical regimes: < m < 7 is the ferromagnetic critical transition, and 7 < m < tt is the 'non-physical 
self-dual line', governing the critical Berker-KadanofF phase |32j . 

In the present paper, we are interested in the case a = 7r/2, which we call the Z2 staggered model. The parameters 
t!i,D2 are then given by: 

vi _ sinu V2 _ cos(7 - u) (2 11) 



^ sin(7 - u) ' 
There are again two regimes: 

• Regime I: 7 < u < 7r/2. 

In this regime, we have ui < and V2 < 0. This was identified by Baxter [5] as the location of the antiferro- 
magnetic critical transition. There is an isotropic point u = ^ + ^ for which vi = V2- The exact solution on the 
lattice, and the corresponding field theory were studied in detail previously [31 13]. 

• Regime II: < u < 7. 

In this regime, we have now wi > and V2 < 0. There is thus no isotropic point, but for u = ^ we have that 
vi = —V2 differ only by a sign; we shall see that at this point isotropy is recovered in the continuum limit. In 
the following, we will focus our attention on this regime. 
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Figure 5: Lines of spectral parameters for the staggered model. The block 7?.-matrix is represented in bold lines. 



Let us now recall the lattice structure of the staggered model [3l [4] . This is better described by the block 7?.- matrix 
(see Fig. [5]): 

7^JJ + l(li) = R2j,2j+l{u - 7r/2)^2j-l,2j(M)^2i+l,2j+2(w)^2j,2j + l(w + 7r/2) ■ (2-12) 

Using (|2.9p . we get the expression for TZjj+i{u) in terms of the TL generators ej: 

TZjj+i{u) = — - sin^(27 — 2u) 1 — - sinwsin(27 — 2w) [cos(7 — w)(e2j-i + e2j+i) + 2cos7COSu e2j] 
+ i sin2wsin(27 - 2u) {e2j~ie2j + e2je2j_i + 62^62^+1 + 62^+162^) 

+ sin^ wcosw [cos(7 - u) (e2j_ie2j+ie2j + e2je2j_ie2j+i) - cosu 62^62^-162^+162^] . (2.13) 
As a consequence of the Yang-Baxter equations, 7?.j.j+i(u) commutes with the operator: 

Cj = (cOS7)-2 7?2j-l,2j(V2)i?2j + 1.2j+2(V2) • (2.14) 

Furthermore, using the TL algebraic relations, we see that c| = 1. Thus, the 7?.-matrix has a Z/2Z symmetry, 
arising from the staggered structure. We call t(u) the one-row transfer matrix, with horizontal spectral parameter u 
and vertical spectral parameters (0, 7r/2, 0, tt/2, . . .). We consider a lattice of width 2N sites and height 2M sites. The 
partition function of the staggered model with periodic boundary conditions is then: 

Zmn = Tr [t{u)t{u + TT/2)f' . (2.15) 

The two-row transfer-matrix t{u)t{u + 7r/2) commutes with the Z/2Z 'charge operator': 

N 

C ^ n [('^0S7)"^ i?2,_i,2,(7r/2)] . (2.16) 



2.3 Very anisotropic limit 

An interesting aspect of Yang-Baxter integrable statistical models is that the transfer matrix generally possesses 
a very anisotropic limit, where its derivatives with respect to u are local, one-dimensional Hamiltonians. In the 
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expression (|2.13p for the block 7?.-niatrix, we observe that T^-jj+i (0) = — isin^27 1, and so the two-row transfer- 
matrix reduces to a cychc translation of two sites to the right in the limit u — > 0. The first-order Hamiltonian 
is: 

N 



2 du 
which gives; 



isin27 ^7e,,+i(0)-i^^(0), (2.17) 



2N 

i7 = 2iVcos27 1 + ^ (-2 cos 7 + ejCj+i + Cj+iCj) . (2.18) 
i=i 

For a generic value of 7, the TL algebra can be represented as acting on a chain of spin-i variables with Ug(sl2) 
symmetry (for q ~ Indeed, consider the Hilbcrt space [span(| t), | i))]*^^^, describing 2N spins (cri, . . . ,<J2n), 

and take two consecutive spins cTj, cTj+i'. the total spin <Tj -I- cTj+i can be in the representation of spin or 1. We call 

Pj^j^i the projector onto spin-0 according to this decomposition. Then it turns out that the unnormaliscd projectors 

ej- = 2cos7 Pjj+i (2.19) 
satisfy the TL algebra (|2.6p . In terms of the Pauli matrices, the generators (|2.19[) can be written: 

Using this representation, the first summand in (|2.18p would simply give the XXZ spin chain. The quadratic terms in 
Bj, Cj+i are due to the staggering, and lead to a different spin model. After some algebra with the Pauli matrices, we 
obtain: 

-e-"^^^ {<J+^,<y-+2 + '^7+i<+2) - ^^I+i'^I+2 • (2.21) 
Let us write our quadratic TL Hamiltonian ()2.18p more generally as: 

2N 2N 

H{K,,K2) ^K^Y. + ^2 ^(e,ej-+i + e,+ie,) . (2.22) 
i=i i=i 

Then, from ((2?20)) - ((2?2T1) . the above Hamihonian can be written: 

fj 1 
H{Ki,K2) = I Y'^J- ■ + JiHcj+a-^i + '^74+1) + + iMa]_^ - a|+2)(<'^7+i + <^J 4+1^ ' 

(2.23) 

where: 

= - (i^i + 2 cos 7 i^2) , J^ = -(l COS 7 A'l + , J2 = , J3 = K2 sin 7 . (2.24) 
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Figure 6: Zig-zag spin chain with Ji and J2 interactions. 



The J2 term in equation ()2.23p represents two XXX spin chains, hving on the even and odd sites respectively. The 
Ji terms correspond to an XXZ interaction with a 'zig-zag' shape (see Fig. |6]). The J3 term is an anti-Hermitian 
three-spin interaction, with no obvious physical interpretation. 

In the case of our Z2 staggered model, wc have {Ki,K2) = (— 2cos7, 1), and so J2 = 1, Jf^ = 0, Jf = — sin^ 7. 
Thus the Hermitian terms of (|2.23p correspond to two antiferromagnetic XXX spin chains with a ferromagnetic — cr^cr^ 
zig-zag interaction. Let us discuss the consequences of a naive bosonisation [24] argument for this model. We can 
discard the irrelevant J3 term, and we end up with two free bosons (f>i , 02 with compactification radii Ri = R2 = R 
independent from 7, coupled by the quadratic term J^d(j)id4)2- Since this term is symmetric in the exchange of the 
bosons ((f>i ^2), the symmetric and antisymmetric combinations <i> = ((/>! + ^2)/ V^, 4> = {4>i ^ 4>2)/ 'v/S are decoupled 
free bosons, with compactification radii R± = R± 6R, where SR depends on 7 through Jf. Thus, we obtain two 
decoupled free bosons $, (j> with both radii depending on 7. However, using the Betlie Ansatz exact solution of the 
staggered model (see Section[3]), we find that the continuum limit consists of two free bosons with one radius depending 
on 7, and the other radius independent of 7. It seems then that the bosonisation approach misses an important effect 
due to the anti-Hermitian J3 term in (|2.23p . 

Let us go back to the model (j2.22[) with general values of A'l , K2 , which was studied numerically in Ref . [25] . 
Additionally to the Z2 staggered integrable model, this Hamiltonian contains a remarkable point at {Ki,K2) = 
(—4 cos 7, 1). This was identified in Ref. [20] as the g-deformcd version of the Majumdar-Ghosh (MG) spin chain [19]. 
The latter was defined as an SU(2) spin-^ chain with nearest- and next-nearest-ncighbour interactions, for which the 
totally dimerised state (where the pairs of spins a2j-i,o'2j form a singlet) is the exact ground state. The excitations 
of the isotropic MG chain are known to be gapped spinous [26]. In Ref. [25] it was proved that for 7 = 7r/i with i > 5 
integer, the dimerised state is still a ground state, and numerical evidence was given for the existence of a gapped MG 
phase in the model (|2.22p for generic 7 in the range {n* ~ 1.5) < 2 cos 7 < 2. 

2.4 Relation to anyonic chains 

We end this Section by mentioning an additional motivation for the study of TL Hamiltonians such as (|2.22|) . which is 
their connection to 'anyonic chains', an object introduced and studied recently [HK^^]. The starting point is the Read- 
Rezayi construction [27] of trial wavefunctions for the Fractional Quantum Hall Effect. Following this approach, the 
trial wavefunctions are given by the correlation functions of the Z^, parafermionic CFTs. The elementary excitations 
above the ground state (anyons) then obey the fusion rules given by these CFTs. 

The anyonic chain, defined in Ref. [H] for the Z3 parafermionic CFT, consists of M + 2 anyons with short-range 
interaction. Each anyon lives in a Hilbert space whose basis elements are labelled by the primary fields of the CFT. 
In the case of the Z3 parafermionic CFT, the primary fields are 1 and t, with the fusion rule: 

TXT=1 + T. (2.25) 

We fix M + 2 anyons along a chain to be in the non-trivial state r, and give an interaction energy to the different ways 
they can fuse with one another. Suppose the r-anyons are allowed to fuse according to the linear diagram shown in 
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Figure 7: (a) Anyonic chain, (b) Transformation of fusion diagrams through the _F-matrix. Figures reproduced from 
Ref. ini. 



Fig. [7^, where the Xj are also anyons. LocaUy. the fusion diagram between two consecutive r-anyons and the anyons 
Xj-i, Xj+i can be rewritten in terms of a new anyonic variable Xj, using the i^-matrix (see Fig. Wp)- The variable Xj 
now results from the direct fusion of the two consecutive r-anyons. We give energy — e to the configuration Xj = 1, and 
energy to the configuration Xj = r. Let us describe the Hilbert space and the interaction Hamiltonian in terms of the 
intermediary states Xj. The basis of the Hilbert space is labelled by the words xi . . . xm-i on the alphabet {1, t} which 
are allowed by the fusion rules. These are the words which do not have the word 11 as a subword. Fig. \7i) defines a 
local change of basis \xj) — + \xj). In this new basis, the interaction term defined above is proportional to the projector 
Pj = |1)(1|. Transforming back to the Xj basis, we get the operators pj, given by {x'\pj\x) = {Ftx-+1)x'i{Ftx-+1)ix- 
The total Hamiltonian is then defined as 

A/-1 

-Canyons = X! ' (2.26) 

3=0 

We are now ready to state the correspondence between the anyonic chain and our Hamiltonian H (|2.18p . The 
expression of the i^-matrix is such that the unnormalised projectors ej = (2 cos f )Pj satisfy the TL algebra (|2.6p for 
y/Q = 2cos^ [21] • In fact, the allowed words xi ...xm-i are in bijection with the A4 RSOS configurations, and 
the operators ej act on them as in the A4 RSOS model. This exact correspondence holds in general between the 
anyonic chain and the Ak+i RSOS model, with the loop weight y/Q = 2 cos 

Let us recall briefly the construction of the RSOS representation of the TL algebra [28], in the case of the Ap 
models. The basis states for the vector space ?i are labelled by M = 2N height variables {hi, ... , /i2Af), such that, for 
all j: 

h, e{l,...,p}, \hj-h,+i\ = l. (2.27) 
Then we define the operators {ej)j=i,,,2N acting on Ti. by their action on the basis states: 



Jsihjjsih^) 

ej\hi...h2N) = —\h,...h'^...h2N), (2.28) 

S{h) = sin^^. (2.29) 

p+l 

It can be shown [55| that these operators satisfy the TL algebra (|2.6p with loop weight y/Q = 2cos7r/(p +1). In the 
particular case of Z3 anyons, the bijection xj hj reads: 

1 ^ 1 , r ^ 3 for j odd 

1^4, r ^ 2 for j even 
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Note that this maps the anyonic words xi . . . xm-i to the subset of RSOS configurations such that even sites j carry 
even heights hj. However, this subset and its complementary (where even sites carry odd heights) are not coupled by 
the ej. The two sectors are related by a one-site translation. Hence the equivalence between Z3 anyons and A4 RSOS 
is valid up to a degeneracy factor of two. 

From this equivalence between anyonic chains and RSOS models, one understands that -ffanyons (|2.26|) is just an 
RSOS version of the well-studied XXZ model with anisotropy A = \fQl2. It is known that the corresponding 
effective field theories arc the CFT minimal models. More interestingly, in Rcf. [22j, a three-site interaction was 
introduced: 

M-\ M 

-^anyons = " XI + ^ H Pj" ' (^'^^^ 

where p'^ denotes the projector onto the 1-channcl for the fusion of the r-anyons at positions [j, j + l, j + 2). Hamiltonian 
-f^anyons ^hc RSOS representation of the quadratic XL Hamiltonian (|2.22p . In the context of anyons, only Beraha 
values \/Q = 2 cos7r/(p -I- 1) with p = 3, 4, 5, . . . are allowed, whereas the loop formulation exists for generic Q. Hence 
H{Ki,K2) (|2.22p is really a generalisation of -ffanyons (|2.30p . Numerical studies [HI [23] have shown that the phase 
diagram of this model is very rich. In particular, it contains a critical phase rcnormalising to the XXZ Hamiltonian, 
a gapped phase of MG type (see above), and another critical phase governed by the integrable Z2 model which is 
studied in the present paper. It was also shown in these works that the transition between the XXZ and MG phase is 
in the universality class of the dilute 0(n) model, with n = \fQ- We recall the results of |5S] in Fig. [51 where we have 
used the same conventions as in |22| to parameterise the interaction: 

Kx = 2\/^ sin0 - COS0 

The phase diagram is represented in the range [0z2jO]i where = Arctan Q^^''^ — tt is the value corresponding to 
the Z2 model. 



3 Bethe Ansatz solution 

In this section, we present the solution by Bethe Ansatz of the model presented in Section [51 We find that the Bethe 
roots form two coupled Fermi seas, and the elementary excitations are holes close to the Fermi levels. In the continuum 
limit, we obtain the dressed momentum and energy (j3.2ip of the holes, and the dressed scattering amplitudes ()3.26p 
between them. The central charge of the theory is c ~ 2. Using the Wiener- Hopf technique for the computation of 
finite-size corrections (see Appendix B), we derive the conformal spectrum ()3.30p . It has the form of a two-component 
Coulomb gas. 



3.1 Bethe Ansatz Equations 

We use the Algebraic Bethe Ansatz [1], and we define the Bethe roots Uj as: aj = 1(7 — 2uj). The Bethe Ansatz 
Equations (BAE) and the eigenvalues and eigenvectors of the transfer matrix t(u) in the 7'-particle sector are: 

^ ^ sinh \ {aj — ai — 2i'f) 



sinh(aj — 17) 
sinh(aj -I- 17) 



1=1 



■ ai + 217) 



(3.1) 



e'^[sin2(7-.)]-n ^^";;!i^^+;"-+^^";i 

fj-^ smh ^ [17 - {aj + 2m)] 



a,- 



- 2iM - 217)] 



e-'n- sin2.)- n ""^tf^"; ^ ^ 

sinh i [i7 + [uj + 2m - 217)] 
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n* ~ 1.5 2 

Figure 8: Phase diagram of the quadratic Temperley-Lieb Hamiltonian (|2.22p . using the parameterisation (|2.3ip . The 
shaded area is the gapped phase, governed by the MG hne (dotted hne). Full lines represent phase transitions from 
the MG phase to the XXZ critical phase (bottom) and the Z2 critical phase (top). The Z2 staggered model sits on 
the topmost horizontal line 9 = ■ 



\-^{ui,...,Ur))=B{ui)...B{Ur)\0). (3.2) 

In Eq. (|3.2p . wc have used the notations of JJ for the monodromy matrix elements. The Bcthe states |5'(ui, . . . , Ur)) 
are invariant under the two-site cyclic translation c~^'^. In the very anisotropic limit u ^ 0, the transfer matrix 
becomes: 

^(0)t(|) = (-^)"e--, (3.3) 
and the corresponding eigenvalue is, from (|3.2p : 

A(0) A f ^) = e^'^ f-'-^) ff ^^^Ma^+^^l . (3.4) 

The energy for the Hamiltonian (|2.18p is the logarithmic derivative of the eigenvalue: 

1 . dlog[A(u)A(u + 7r/2)] ^ 2sin2 27 



E= — sin 27 



2 du 



= 2N cos 27 - y , . (3.5) 

u=o ^ cosh2aj -cos27 



Equations (|3.4p - (|3.5p show that each Bethe root a contributes to the total momentum and energy by: 

o, -1 sinh(aj-i7) 2sin2 27 

2fc, = — ilog . , , — , e, = — — . 3.6) 

smh(Q;j + 17) cosh 2aj — cos 27 

Because of the periodicity property of the Boltzmann weights R{u + tt) = —R{u), the Bethe states are unchanged 
under a — > a + 2i7r for any of the Bcthe roots. So we can restrict our study to the strip: < Im a < lix . The root a 
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gives a negative contribution to the energy (|3.5p if it is of the form: 

=Af or a«=A«+i7r, (3.7) 
with A^"'' reaL So, at low energies, the system is described by two coupled Fermi seas 

{Xf^A^f^}- The BAE for the 

A^ ' are: 

2Nk{\f^) = 2TTlf + 20 - ^ ^ e(''-^)(A^"' - \f^) , a = 0, 1 , (3.8) 

6=0,1 1=1 

where r'°) is the number of roots a^"-* . The momentum, energy and scattering phases arc given by: 



(3.9) 



„,/>x •-, sinh(i7-A) ^rnWxN • -, sinh(i7 + ^) 

smh(i7 + A) smh(i7 - ^) 

e(A) = -sin27 x 2fc'(A) = , e(±^)(A) = -ilog ^"^^^^^ + f j . 

^ ' ' ^ ' cosh2A-cos27' ^ ' ^ cosh(i7 - |) 

The Bethe integers satisfy /j""* € ^(iV + r*^") — 1) + Z. The total momentum and energy are: 

= |EE(^r + ^)+-(-^"^+^^^^)' (3-10) 

a=0,l j=0 ^ ^ 

^ = 27V cos 27+ ^ E^^^^)- (3.11) 

a=0,l j=l 

A special case is when the Bethe roots {A^"^}, {Aj-^"*} are identical on the two lines: we call these symmetric states. 
It is a remarkable fact that this subset of the spectrum is exactly the complete spectrum of the XXZ spin chain on a 
periodic lattice with N sites: 



ffxxz = -IJ2 + «+i + ^0 «„+i] , (3.12) 

?n — 1 

with Ao — — cos 27. Indeed, we have the identities: 

2fc(A) ==fcxxz(A), e(A) = exxz(A), (6^ + e(±i))(A) - exxz(A) , (3.13) 



where quantities with the subscript 'XXZ' are related to the Bethe Ansatz for XXZ. Thus, the BAE (|3.8p for symmetric 
states are equivalent to the XXZ ones, and the energies are related by = 2i?xxz- 

3.2 Continuum limit 

The ground state corresponds to r(°^ = r*^^^ = N/2, with the Bethe integer distribution (see Fig. [9^): 



^(0) ^ ^(1) ^ N/2+1 



2 +j, j = l,...,iV/2. (3.14) 
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The continuum limit is defined as: 

N^oo, r'^"^ /N -*l/2. (3.15) 

In this hmit, we assume that the spacing between Bethe roots scales like 1/N, and we describe the Bethe root 
distribution by the densities: iVp(°)(A^"^) = 1/(A^"\ - A^"^). We denote [-C^°-\ B^"-'^] the interval spanned by the 
roots A^"''. The BAE equations p.Sp become Lieb equations for the root densities: 

2fc'(A) = 27rp('^)(A) - y / d^i p'^^Xi^L) K'^'^-^XX ^ ^i) , a = 0, 1 , (3.16) 

where the kernels are given by: = (e^'^))'. In the ground state, we have C^^'^\B'^^'^^ oo, so Eq. (|3.16p can 

be solved by Fourier transform. The solution involves the symmetric and antisymmetric inverse kernels J*-*-* (see 
Appendix A): 

1 + J(±)(c^) = ^ — ^ . (3.17) 

2ti -[Ki°){uo)±K(^){uj)\ 

The ground-state densities are p*-"-* = p^^-* = p^o, where: 

Poo = (<5 + j(+))*(2fc')/(2^)-T ,y (3.18) 

47 cos[7rA/(27)J 

The symbol 7k- denotes convolution. 

An elementary excitation above the ground state consists of a hole Xh in the distribution {x'f^} or {A^^'}, interacting 
with all the particles in both Fermi seas. Let A be a physical quantity defined as: 

^ = ^ E E"(^5"^)^ / dXp^'\X)a{\)+ [ dApW(A)«(A). (3.19) 

In the presence of a hole X^, the variation of A with respect to the ground-state value ^0 is given by the dressed 
quantity ad (see Appendix A): 

A - Ao = ^ aaiXh) , = -{S + J^+^) * a . (3.20) 



The momentum and energy of a hole are thus: 

2fcd(A) = -2Atan 



tanh ( fA 



(3.21) 

^'^^ ' 27Cosh[7rA/(27)] ' 

In the region A — > 00, the dressed momentum is close to the value — 7r/2, and the dispersion relation is linear, with 
Fermi velocity v: 

TT \ TT sin 27 



tdc^-v {2kd + -j , v= ' . (3.22) 



Hence, hole excitations are gapless, and the theory is critical. 
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3.3 Dressed scattering amplitudes 

In the presence of holes, the root densities p^"^ coexist with the densities of holes p^yf' . The Lieb equations (|3.16p 
become: 

\ 2^(p(i) = 2fc' + /v(i)*p(")+i^(o)*pW. 
These coupled equations can be rewritten in terms of p^""^ + jo|°^ and 



(3.23) 



(3.24) 



where the kernels 7*^°^ J'^^-'^^ are defined as: 



j(°)^i(j(+)+j(-)), J(±i)^l(j(+)- J(-)). (3.25) 



The Fourier transforms of the kernels J^^-', J^^^^ are: 

sinh(7r — 37)0 



J(0)(a;) 



J(±i)(u;) = - 



2 cosh7w sinh(7r — 27)0; ' 
sinh 7C1J 



(3.26) 



2 cosh7w sinh(7r — 27)0; 



3.4 Central charge and conformal dimensions 

The low-energy spectrum consists of 'electromagnetic' excitations above the ground-state distribution, similarly to 
the XXZ case |30| . In the present case, the Bethe integers distibutions {lf^},{lf^} 

can be chosen independently, as 

depicted in Fig. [9)D-[9t. A magnetic excitation consists in removing m'°) (resp. m^^^) roots of type (resp. ot"^) 

from the ground state, while keeping the Bethe integer distrubutions {if*'^^} symmetric around zero. An electric 

excitation consists in shifting all integers if^ (resp. by e'"' (resp. e'^'). 

The central charge and critical exponents are obtained from finite-size corrections to the total energy and momen- 
tum: 

TTC 

Eq ~ Neoo - " X ^ . (3-27) 

E-E, . .xM^, (3.28) 
27r(A-A) 

Q = (3.29) 

where Eq is the ground-state energy. 

We start by discussing the untwisted case </) = 0. The ground state is a symmetric state, and thus the ground-state 
energy is twice that of the XXZ spin-chain (|3.12p . The latter has a Fermi velocity v and central charge one. Using 
Eq. p.27p . the central charge of the staggered model is then: c = 2. The finite-size corrections to the energies for the 
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X X X X X X 



-3-2-10 1 2 3 



Ij ■ X X X X X X 

(a) 



r(l) 

Ij ■ XXXXXX XXXXXX 



-3-2-10123 -3-2-10123 



r((>) 

1] ■ XXXX XXXX 

(b) (c) 

Figure 9: Example of Bethc integer distributions for = 12. (a) Ground state; (b) Magnetic excitation m'"' = 2; (c) 
Combined magnetic-electric excitation trS^^ — 2,e'"-' = 1. 



electromagnetic excitations are computed from the Bethe-Ansatz solution in Appendix B. They yield the conformal 
dimensions: 

1 / e , — \ 1 



Aem,em = -( ^= + 777.^2.9 j + -{e + mY 



1 / \ 2 1 

1 / e ; — \ 1 



Ae™,e™ = -(^-mv/2gj +-(e-m)^ (3-30) 



TT - 27 ^ 1 

' - «<^<2' 

where: 

(3.31) 

e = e^"^ — e'"'"'' 777 = 777'*'-' — 777^^-' . 

When the twist (j) is not zero, the above exponents are still correct, with the change e — > e + 2(/)/7r. In particular, 
the staggered Potts model corresponds to a twist (jj = j = ttcq. The groimd state has an electric charge e = 2eo, with 
exponents A = A = eQ/(4g), so the effective central charge is: 

ctw = 2 , Co = - = --5. (3.32) 

9 TT 2 

3.5 Application to the calculation of critical exponents 

In the loop formulation, we obtain the fc-leg dimensions as follows. For any system size N, the number of legs k must 
be even, and the conformal dimensions are defined with respect to the twisted ground state: 

hk^hk = Ak-^. (3.33) 
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The fc-leg dimension corresponds to a magnetic defect m = fc/2, with a minimal value for m and electric charges 
e = e = (no background charge). There are then two distinct cases: 

^ - ^ if fc ^ [4] 

^2 2 (3.34) 

^ + iffc^2[4]. 
16 8 45 ^ ^ 

Similarly, the magnetic exponent of the staggered Potts model is defined with respect to the twisted ground state: 

hH^fiH^AH-^. (3.35) 
4.9 

The magnetic dimension Ah corresponds to a twist = 7r/2, which forbids any non-contractible loop around the 
cylinder. Before we obtain A^, we need to discuss the conformal dimension for the sector m = = with a general 
twist (f). In the regime < < tt, the lowest dimensions are: 

Ai((/)) = A(0/^_o),(o,o) = -^^^^ , (3.36) 

A A (0A - 1/2)^ 1 

A2(0) = A(0/^_i^o),(i,o) = 1^ (3-37) 

A3(0) = A(0/^„2,o),(o,o) = -^^-^^^ ■ (3-38) 

The lowest dimension is respectively Ai, A2, A3 on the intervals [0, (po], [<f>o, l — ipo], [1 — (^0, i"], where ^0 = ''■(l + 2,9)/4. 
See Fig. [TOl In particular, for = tt/2, we get Ah = A2(7r/2) = 1/8, and thus: 



3.6 Numerical checks 

We have verified the above expressions for the effective central charge ctw and the fc-leg exponents hk by numerical 
diagonalisation of the transfer matrix at the pseudo-isotropic point u = ^ . 

As usual, the critical exponents can be extracted from the finite-size corrections in TV to the dominant eigenvalues 
in the various sectors labelled by k. We consider the geometry of a strip of width 2N strands with periodic boundary 
conditions in the transverse direction. Estimates for ctw (resp. hk) are then obtained from fits involving three (resp. 
two) different sizes N. We use even N throughout. Odd N introduces a twist that leads to different effective exponents 
that we do not consider any further. 

It is convenient to parameterise 7 = j through a new parameter t. The results for ctw with iV < 14 are shown 
in Fig. [11] The agreement with (j3.32p is excellent. Results for Xk = 2hk with < 12 are given in Fig. [12] for 
k = 2,4:, 6, 8. The agreement with p.34p is very satisfactory, especially when fc = [4]. 



4 Toroidal partition functions 

In this section, we use the results from Section [3] to construct explicitly the continuum partition fimction Z of the 
statistical model on a torus. Assuming that the conformal spectrum p.30|) obtained from the analysis of the BAE 



17 




7r/4 7r/2 37r/4 tt 

Figure 10: Lowest conformal dimensions in the twisted sector, for 5 = 1/4. 




1 1 I \ I \ I \ I \ I I 

'3 4 5 6 7 8 

t 



Figure 11: Numerical estimates of the effective central charge ctw, as compared with the exact expression p.32p . 
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t 



Figure 12: Numerical estimates of the fc-leg exponents X = 2hk with k = 2,4,6,8, as compared with the exact 
expression p.34p . 



is complete, we sum the conformal characters over all possible conformal dimensions to obtain Z . The resulting 
expression (|4.14p for Z shows that the continuum limit of the model consists in one boson and two Majorana fermions, 
which decouple in the bulk and couple only through boundary conditions. We discuss only the untwisted case here, 
leaving the twisted and Potts model cases (including the study of particular values of Q) to Appendix D. 

We denote by r the modular ratio of the torus, and we write q = e^"^"^. Other notations are defined in Appendix C. 
The primary states of the corresponding CFT have conformal weights Aem,?m and Aem,?m given by the Bethe Ansatz 
results p.30|) . where the charges satisfy the parity conditions: 

e + ee2Z, m + 7ne2Z. (4.1) 

The partition function on the torus is given by the sum of the generic conformal characters xa,Xa- 

Z{g) = Tr ^-£o-c/24^ ^ ^ XA(g)XA(9) , (4.2) 

A, A 

where the sum is over all possible primary states, and xa is the trace of (j^o^'=/24 over the descendants of the primary 
state $a: 

XA(g) = TrA <z^«-^/^^ (4.3) 

The character xa can be inferred from the possible Bethe integer distributions. Starting from an electromagnetic 
excitation with dimension A, we can create vacancies, by shifting the largest Bethe integer Ij Ij + n,n > 0. This 
vacancy state has dimension A + n. These vacancies can be combined, and the state with shifts (rii, . . . ,71^) has 
dimension A + ni + ■ ■ ■ + rik- Furthermore, vacancies can be introduced independently on the two lines Ij^^- Let 
us denote p{n) the number of partitions of the integer n. We have: 

XA(g) = E p(n(°^Mn^^^)9^+"*"'+"'''-^/^^ = ^^^S^ > (4-4) 
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where 77(t) is the Dedekmd function (|C.3|) . Using (|4.4p with c = 2 and the parity conditions (|4.ip . we obtain: 

/ \ 

1 



Z{9) 



Hr)V 



E 



E 



E + E 



m,m even m,m even m,m odd ra^m odd 

e,e even e,e odd e,e even e,e odd y 



(4.5) 



Using the Poisson summation (jC.Sp . this can be written: 



m,m' even 



where: 



A 



E 



ra,m' even 

^ Z„,^,„,(l/2) 




(4.6) 



(4.7) 



C — ^ ^ -Z^m.jTi' (1/2) , D — ^ ^ Zjn^m' 

rn odd,?n' even m,m' odd 

and Zm.m'ig) is the bosonic partition function with defects m,m' (see (|C.4[) '). 



The partition sums A, B, C, D can, in turn, be expressed in terms of the Jacobi ones (jC.9[) . using (jC.Sp again: 
A = \{Zl + Zl + Zl), B = i(-Z| + Z| + Z|), 

D = \{Zl-Zl + Zl). 



c 



zl-zl), 



(4.8) 



Using the transformation of Jacobi and Coulombic partition functions under modular transformations, one can show 
easily that the expression ()4.6|) is modular invariant. Let Z(r,r') be the partition function of the Ising model on a 
torus with respective boundary conditions on the spins a in the two directions of the torus: 

a^(-l)V, a^i-iy'a, r,r' e{0,l}^ (4.9) 

Using the relation (|C.10[) between Z(r, r') and Z^, the partition sums A,B,C,D are written in terms of the Z{r, r'): 



A 
B 
C 
D 



i [Z(0, 0)2 + Z(0, If + Z(l, 0)2 + Z(l, If] , 
i[Z(0,0)Z(0,l)-Z(l,0)Z(l,l)] , 
i[Z(0,0)Z(l,0)-Z(0,l)Z(l,l)] , 
i[Z(0,0)Z(l,l)-Z(0,l)Z(l,0)] . 



Hence, from (|4.6p and (|4.10l - [4T3|l . the partition function Z{g) reads: 



(4.10) 
(4.11) 
(4.12) 
(4.13) 

(4.14) 



m=ri+r2 [2] 
m'=r[+r'2 [2] 
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The degrees of freedom contained in Z{g) are a compact boson Lp (see (|C.ip ) with couphng constant g = {n — 2^) / {2tt) , 
and two sets of Ising spins ai,a2- The boundary defects for 95, (Ti,(T2 are respectively (m, m'), (ri, rj^), (r2, rj), and 
obey parity conditions, as shown in (|4.14p . Apart from these conditions, the three degrees of freedom (p,(Ji,(j2 are 
decoupled. These results are very similar to what was found in [3T] for a lattice model related to = 1 supcrconformal 
theories, where only one Ising spin was present. 

Like it was done in [3T] for the 19- vertex model, here we can also identify the degrees of freedom (f, cri, (J2 in the 
lattice model. For this purpose, we consider the vertex model defined by the block 7t-matrix (see Fig. [5]). It was 
shown in [5] that there are 38 possible vertices. Each edge can be in one of four states: 1,1, \, \\. Let A^Q(r) be the 
number of edges adjacent to the site r, which are in the state a. An essential property of the model, arising from the 
combination of the magnetisation conservation and Z/2Z symmetry, is that A^^i and N^^ are both even for every vertex. 
Thus, for a given lattice configuration, the lines formed by the | and || edges can be viewed as the domain walls of two 
distinct Ising models, both living on the dual lattice. The remaining edges carry arrows, which define a height (SOS) 
model on the dual lattice. Although these three degrees of freedom are coupled in the lattice model, our results on the 
continuum partition function show that they decouple in the continuum limit, except for their boundary conditions, 
which keep track of the parity of domain walls and arrows around each direction of the torus. 

5 Integrable massive deformation 

In this Section, wc follow the approach of jl3j to construct a massive deformation of the lattice model, and study 

its excitation spectrum. Using the dressed scattering amplitudes, we obtain partly the ^-matrix for elementary 

excitations. Then we use known results from [14j on this S-matrix to conjecture a TBA diagram, and we use the 

TBA equations to calculate the ground-state energy scaling function. In the UV limit, we retrieve the results from 

Section [3l Finally, we use these results to propose an effective QFT for the massive deformation, which is a complex 
(2) 

version of the ' Toda theory 

5.1 Massive integrable deformation on the lattice 

We now consider a deformation of our model where the spectral parameters acquire an extra staggering, this time 
in the imaginary direction. We choose the pattern u + iA/2, u — iA/2, u + tt/2 + iA/2, u + tt/2 — iA/2. This kind of 
construction has been widely used to induce an integrable massive deformation from integrable lattice models |13[ I33j . 
We obtain a modified set of Bethe equations: 

2Tr{p + p/Jf") (A) = 2k' {X + A) + 2k' {X - A) + {K<-°^ ★ P<"))(A) + (A^^-^) * p''')(A) 

(5.1) 

2t:{p + p/Jd) (A) ^ 2fc'(A + A) + 2fc'(A - A) + (if (1) * pW)(A) + (A'W * p(i')(A) . 

To explore the corresponding physics, we write what is often called the physical equations, that is, the equations 
describing scattering of dressed excitations. Since the ground state is obtained by filhng up the p^"^ and p^^^ fines, 
this is easily done by reexpressing the equations so that the densities of holes appear on the right-hand side (see 
Section [331). Wc find: 

2tt{p + p^)(") = s + $("^") * p|f ' + * 

(5.2) 

2^(p + p^)W = s + $(1^0) ^ p(0) ^ ^ ^(1) ^ 
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where = -27rj(''-^), and: 

1 



27 



cosh 5^ (A - A) cosh 5^ ( A + A) 



(5.3) 



This function has tails at |A| ^ A where s(A) decays exponentially as in the masslcss case. These describe a 'ghost' 
of the initial massless theory, whose physics does not depend on A, and which we will not discuss in the following. It 
decouples entirely from the region where |A| <^ A, which is of interest to us. In this region, we have 

s(A)«^e--^/(27)cosh^, (5.4) 
7 27 

so the corresponding momentum and energy are 

2fcd(A) = -/ s(/i)d/^ « -4e-"^/(27) sinh — , (5.5) 
Jq 27 

€d{X) = sin 27 s(A) w4v 0-^^/(2'') cosh— . (5.6) 

27 

We thus obtain a massive relativistic spectrum, as happens systematically in this kind of construction. The mass is 
given by: 

ttA 



^^4exp(^--j. (5.7) 

The question is then, what kind of scattering theory do we obtain, and what quantum field theory does it correspond 
to? 



5.2 Scattering theory 

To answer the above question, we start by reinterpreting the kernels as derivatives of scattering phases between 

basic particles. We will now denote the holes in the 0, 1 sea by the labels 0— and 1—. If we rescale the parameters 
A, uj to: 

TT 



27 



X A . 



I. 27 

TT 



(5.8) 



and set t = n/^, we obtain, up to a constant phase which will be obtained below: 



dk _ sinh^(t-3) 

— sm k& ; 

k sinh^(t- 2)cosh^ 



(5.9) 



S^_°!!>{e) = Sl'!^'ie)o,i exp 



(1,0), 



dk 



■ sin k9 



sinh ^ 



sinh^(t - 2) cosh ^ 



These two 5- matrix elements can be interpreted in terms of the scattering matrix 5y(/3sGi 9) of the Sine-Gordon (SG) 
model [34], with action: 

2' 



-dt^ipd^ip + cos(/3sG'p) 



(5.10) 
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If we set: 



then S^^f^^ is the kink-kink (or antikink-antikink) scattering element for the SG model 



SG 



t - 1 



(5.11) 



(5.12) 



It is natural to identify the holes as two types of antikinks (0, 1). We expect the scattering theory to contain also 
a corresponding doublet of kinks, with a full scattering within the (0,0) and (1, 1) sectors described by two copies of 
the Sine-Gordon S'-matrix. 

We now observe that the kernels s'^!}\ 5*1^1°^ are related to the SG scattering matrix S with an imaginary shift 
in the rapidity 6 [HI. The scattering theory defined by S*'"-") S'^'^' 5, S"'"'^' = S't^^"' = § was introduced in 
Ref. [14], where it was proposed as the scattering theory for left/right (L/R) massless particles describing the flow 
between minimal models of CFT under a perturbation by the $13 primary operator. In Ref. [14j . using the unitarity 
and crossing conditions, the normalising factors for the S-matrix were computed. The resulting scattering theory is: 



Four basic particles 
S(°'i)(^) 

^(1,0) 



0+,0-,l+,l- 

m 



-m 

Z(9-m'-^' 



s(psg;0 



_t-2 
2 

t-2 



(5.13) 



where the normalisation factors read: 



ZiO) 



m 



sinhi^L^ 



coshf^^ 



exp 



■ exp 



dk sinhf^(t-3) 

— sm kO ; 

k sinh^f (t- 2)cosh^ 



(5.14) 



dk 



■ smkO 



sinh ■ 



sinh4?(t 



2 2) cosh ^ 



From (|5.1ip . we see that the SG 5-matrices are in the attractive regime for t £ [2,3] and repulsive regime otherwise. 
We stress that 0, 1 are not antiparticles of each other. 



5.3 Ground-state energy 

The scaling function for the ground-state energy is the relevant object to describe the RG flow of a scattering theory. 
We consider the system on a finite circle of circumference R. Then the ground-state energy E{^,, R) has the scaling 
form: 

Eifi,R)^^F{fiR), (5.15) 



where fi is the mass of the elementary particles, given in (|5.7p . 

In the present case, the ground-state energy can be obtained simply, using the following identity on the dressed 
kernels: 

-(0 1),' >! 27r sinh(7r — 47)ti;/2 



$("^")(^) + '|("'i'(^) 



2 cosh7Cj sinh(7r — 2^)uj/2 ' 



(5.16) 
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The right-hand side of (|5.16p is exactly, in terms of the same rapidity 6, the Sine-Gordon kernel but for yet another 
value of the coupling, given by 



(5.17) 



Stt t 
In other words, 

S°±{0) = 5__(/3sg;(?) . (5.18) 

Assuming that the symmetry is not broken between the two types of roots in the ground state, it follows immediately 
that the ground-state energy (calculated, e.g., by the method of Ref. [33]) is twice the ground-state energy of the 
Sine-Gordon model with the same mass for the kinks, and at this renormalised value of the coupling: 

E{^l,R)^2£{0sG;^^,R)■ (5.19) 

This result is in fact quite obvious if we recall that symmetric solutions to the Bethe equations satisfy precisely the 
same system as in the XXZ chain, whose staggering produces the Sine-Gordon theory in the continuum limitQ 

Of course, it would be more satisfactory to establish the result (|5.19[) directly from the scattering theory. We 
obtain this in the RSOS version of the model (for t integer). The ground-state energy E{^, R) is generally obtained by 
the Thermodynamic Bethe Ansatz (TBA) for relativistic scattering theories, introduced in Ref. [35]. The idea of the 
method is to consider a Euclidean theory on a semi-infinite cylinder of dimensions R x L, and to write the partition 
function in two ways: 

Z{R, L) ^ exp [-~E{n, R)L] = Tr (c"-"^^) , (5.20) 

where is the Hamiltonian on an infinite domain when i — > oo. The problem of computing £'(/i, R) is thus equivalent 
to the computation of the free energy on an infinite domain, at finite temperature 1/ R. So one has to find the density 
of elementary particles which satisfies the B AE (|5.2|) , and maximises the free energy at temperature \/R. This results 
in the non-linear integral equations and the ground-state energy, given in terms of the pseudo-energies ta |35| : 



b 



i. = log(l + e-), 0a. = (5.21) 



E{fi,R) = -Y,Y- I La{9) coshe de, 



where Ha is the mass of particles of type a, and Nab is the adjacency matrix of a diagram describing the scattering 
between particles. In a diagonal (non-reflecting) scattering theory, the /Xa and (pab would be given directly from the 
dispersion relations and the 5'-matrix for elementary excitations. However, the present model does allow reflection of 
the particles. The main difficulty here is then to find the correct TBA diagram and masses for the S'-matrix we want 
to study. Following the ideas of Refs. [TnidlKIl] (see the Introduction), we conjecture that the TBA diagram for the 
scattering theory (|5.13p with mass ji (|5.7p is the diagram of Fig. [Ql 

Now, assuming the above conjecture is correct, we check that the TBA equations (|5.2ip for the diagram of Figure[T51 
lead to the central charge (|3.32|) in the UV limit R^ Q. As in Ref. [15]; the ground-state energy in the UV limit is 



^The careful reader might worry about the role of A in both points of view. The staggering in the equivalent XXZ system involves 2A, 
but the anisotropy is also doubled, so the physical mass remains the same. 
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1 2 3 



t-3 



Figure 13: The TBA diagram for the RSOS restrictions of our theory. The leftmost and rightmost nodes are both 
massive. 



obtained in terms of the hmiting values of in the UV and IR limit: 

"t-3 / \ t-3 



1 
nR 



Xa — lim 



-ea(0) 



Va = lim 



-ea(0) 



where C is the Rogers dilogarithm: 



C{x) = -l I dt 



logi , log(l-t) 



1 - t 



t 



(5.22) 



(5.23) 



The quantities Xa,ya are determined by the adjacency matrix Nat and the masses fj,a |15j : 

t-3 



xl = l[il + x,r-^, yl= \{{l + v,r-^. (5.24) 

fc=l 61^6=0 

To connect this with known results [15j on the RSOS central charge, we introduce the quantities Za which satisfy: 

t-3 

^' = 11(1 + ^'')'^"' « = 2,...,t-3, (5.25) 



and write (|5.22p as: 



i^(.,i?)--^ 



t-3 



a=l 



b=2 



t-3 



a=2 



1 



■t-3 



a=2 



t-4 



a=2 



Va 



Va 



(5.26) 



The above expression is exactly the sum of ground-state energies for the At-2 and At_3 RSOS models, so the central 
charge is: 



t{t-l) 



t(t-2)' 



{t-l){t-2)_ 

which is the central charge (|3.32p of the critical theory. 

Finally, we show that, throughout the scaling regime, the ground-state energy E{fj,, R) is twice that of the 
corresponding twisted Sine-Gordon model. Different cases arise, and we will discuss only one: the case when 
t — 3 = 2n -f 1 , n S N. We can then relabel the nodes on the diagram of Figure [131 so that the n leftmost ones 
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1 2 3 

Figure 14: The TBA for the Sine-Gordon model at coupling Psc/i^^) = + + 2). 

are called 1, . . . ,n, the n rightmost ones fi, . . . , 1, and the middle one 0. The TBA equations ()5.2ip then read: 

Ea = dalfiRC0Sh9 -^^Nab (fl* Lb , O = 1, . . . , 7Z - 1 

b 

to. = Sal fJ.R cosh 9 — ''^^ Nab (fi* Li, a=l,...,7i— 1 

b 

e„ = -(j)-k {Ln-1 + Lo) 
= -cl)* {Lj^+ Lo) 

eo = -0*(i„ + L,-O, (5.28) 
where we have set (t>{6) = l/(27rcosh0). We consider symmetric solutions under the exchange of a and a: 

ea = SainRcosh9 - Nab(f>-k Lb , a=l,...,n-l 

b 

e„ = -0*L„_i-0*log(l+ie"^"/2j _^^lQg|^^_j^-eo/2^ 
y = -(p*Ln 

e-a = ea- (5.29) 
The ground state energy is meanwhile: 

£;(Ai,i?) = -2 Lai9) cosh 9 d9 . (5.30) 

•"^ — ' 27r / 

a=l 

We thus see that our system has twice the ground-state energy of a TBA whose diagram is as in Figure [TH and which 
involves a fugacity for the two end nodes of the fork equal to ±i. This is exactly the TBA for the twisted Sine-Gordon 
model, following the lines of Ref. [3S]. For this value of the twist in particular, the results in Ref. [3S| give the central 
charge (using eq.(21) of [3S], with t + 1 = + 2 the total number of nodes; Aj- = — At-i = i in eq.(20) of [5S]): 

c^l-- ^ — . (5.31) 

(n+l)(n + 2) ^ ^ 

Comparing (|5.27p and (|5.3ip for i — 3 = 2?! + 1, we see that c = 2c. 
5.4 The field theory 

We now try to identify the field theory described by our TBA. This of course involves a bit of guesswork. 
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First, recall we have found that the physical mass of the theory scales as fj, oc e~*^/^ where 7 = n/t. On the other 
hand, we can in general expect that we are facing the perturbation of a model of central charge c = 2 by some operator 
of conformal dimension /i, with action: 



A = AcFT 



(5.32) 



where Acft is the action for the critical UV limit. The dimension of the coupling constant uj is [uj] = R^^^^, and thus 
the mass in the TBA scales as /i oc uj^/i^^^'^\ A detailed look at the microscopic Hamiltonian shows that the bare 
coupling is proportional to e~^. If follows that 2 — 2/i = 2/t, and thus: 



h 



t - 1 
t 



(5.33) 



In general, the TBA approach will give for the ground-state energy £'(/i, R) a scries in /i", where the exponent a 
is given by a = 2(1 — h) if correlations of $/i are non zero {e.g. perturbation by in minimal models), a = 4(1 — h) 
if only even correlators are non zero {e.g. the Sine-Gordon model). Other possibilities exist, e.g. a = 8(1 — h) if only 
correlators involving a number of operators multiple of four arc non zero. 

Since the ground-state energy of our model is twice that of the SG model at 0^q/{^i^) ~ (t — 2)/t. it follows that 
a = 4(1 — /).sg)j where /isG = {t ~ 2)/t. Meanwhile, we have identified earlier the dimension of the perturbation as 
h = {t — l)/t, and thus a = 8(1 — h). We arc forced to conclude therefore that in our problem indeed, only correlators 
involving a number of operators multiple of four are non zero. 

Meanwhile, the structure of the scattering matrix suggests the same quantum group symmetry as the one in the 
/3sG theory, with, for generic values of t, only one conserved charge, since the S'^^ elements allow reflection of charges 
between and 1 sectors (this docs not occur at the special points where </(< — 2) is an integer). Finally, the structure 
of finite-size effects showed that the CFT was made of a Dirac fermion and a boson of i-dcpendcnt radius. This all 
leads us to propose that the action of the theory is: 



A[(^,l/'l,l/'2] 



■ 



2) 



(5.34) 



where ip is the boson, ipi , ^2 are the two Majorana components of the Dirac theory. One can check that this theory is 
indeed integrable using the non-local conserved charges ^le"^"'"^/'^^'^ and '02e'-"'"^/'^'''^. The algebra satisfied by these 
charges leads to SU(2)q with quantum-group deformation parameter [37j : 



- exp(— i7r/(5) 



2/32 



47r-^2 



Meanwhile, the basic SG 5'-matrix with the foregoing value of /3sg has also quantum group symmetry 
deformation parameter that corresponds to: 



SG 



87r-/32 



= t- 2. 



(5.35) 
, with 

(5.36) 



SG 



By requiring that the symmetry of the S'-matrix is the symmetry of the action, we identify the two above expressions 
for (5, and we get: 

P"^ t-2 , ^ 

^ (5.37) 



87r 



2t 
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The dimension of the perturbation is indeed: 



1 t - 2 i - 1 , , 

/i = - + — — = , 5.38) 



and clearly only correlators involving a number of operators multiple of four are non zero. We also obtain a non-unitary 
theory, which is expected from the presence of complex terms in the Hamiltonian. 

An important check of our proposal would be to see if the ground-state energy of theory (|5.34p is twice the ground- 
state energy of the related Sine-Gordon model. One might first tackle this question in perturbation theory. We will 
leave this for future work, and content ourselves by examining the question in the limit /? 0. Then /3sg ^ and we 
expect, on the one hand, the ground-state energy to be twice the one of a free boson. On the other hand, our action 
reduces naively to two identical massive Majorana fermions. In the limit /3 ^ however, counter-terms are needed 
and a ip'^ term also appears (exactly as in the case of = 1 theories [MDj leading to an additional free massive boson. 
Denote the ground-state energy of such a boson, and Ef the ground-state energy of a free Majorana fcrmion: 

^'^ J — oo 

/oo 
\og{l + e-^'"''°'^^)cosh0 de . (5.39) 
-oo 

We have the identity 

2Eb{fi, R) = Eb{2^i, R) + 2Ef{fi, R) , (5.40) 

so we see indeed that the ground-state energy of our field theory will be twice the ground-state energy of the SG model 
in the limit of vanishing coupling provided the mass terms are in the proper ratios. More precisely, the left-hand side 
of (|5.40p corresponds to twice the ground-state energy for the theory: 



AM 

so near /3sq = we will need an action of the form: 



I Qa,(^5,^-f Ai^^^^ , (5.41) 



^d„Lpd^Lp + i('(/'l^'01 + V-'2^V'2) 



d^x + Aio J [iJii^ie'^^ + iJ2i^2e-'^^] d^^+^ J cos{2(3^)d^x 



(5.42) 

f 2) 

We now observe that our theory is identical to the Toda theory (more precisely, we need in fact to set n = 1 in 
the more general C^+i theory, whose form is valid for n > 1 only), whose Lagrangian would read [IDl I41j : 



Vi, 1^2] = \d,vd,^ + i(^i#i + ^2^^2) - Mo(^iV'ie^'^/^ + ^2^20-^'^/^) - ^ cosh(V2(7<^) . (5.43) 
2 

Clearly we have to set g ~ i\/2/3. We then see that in the limit /3 ^ the boson has mass parameter twice the one of 
the Majorana fermions, in agreement with (j5.40|) . To summarise our results: 

• The continuum limit of our lattice model is the complex c'^^ theory (|5.42p . 

• The ^-matrix of this theory is given by ((57TO)) . with PIq/{^it) = 2/3V(/3^ + 47r). 

• The ground-state energy is twice the ground-state energy of the Sine-Gordon theory at coupling /3gQ = 20^ . 
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6 Conclusion 



From the point of view of integrable statistical models, one can think of several ways to generalise the construction of 
the Z2 model. First, one can build a model with a staggering of period n > 2, which has a TLjnl symmetry [42]. We 
guess that the corresponding continuum limit will be related to a product of n copies of SU(2)i, made anisotropic by 
a term as in the case n = 2. What the integrable massive deformation might be is however more mysterious. 

Also, the effective field theory for the analog of the non-compact regime [3] is less clear. Another interesting direction 
is to apply the same kind of construction to other models than the six-vertex model. Of particular interest here would 
be the 'dilute version', obtained by staggering the Izergin-Korepin 19-vertex model. 

In the CFT perspective, the expression for the toroidal partition function in terms of Coulombic partition functions 
generally leads to a classification of new minimal series of CFTs. It is possible, in principle, to follow this program in 
the case of the Z2 model partition functions. 

There are also some important questions about the physical interpretation of the Z2 Hamiltonian as a zig-zag spin 
chain. We have seen that non-Hermitian terms in the Hamiltonian play a role, but it could be that the model is in the 
same universality class as a well-defined, Hermitian spin-chain model. Additionnally, at the Majumdar-Ghosh point, 
the gapped excitations above the ground state (spinous) could be studied more systematically, through a variational 
approach similar to [26| . 

In the context of anyonic chains, the RSOS version of the Z2 model is an integrable point in the phase diagram 
of the three-anyon interaction Hamiltonian (|2.30[) . It actually governs the behaviour of this system throughout a 
whole critical phase, as was shown numerically in Ref. |25| . Various features of this phase diagram still lead to open 
questions, such as the complete RG fiow of ()2.30p and the associated operators at the fixed points, but also the 
differences between the RSOS and loop formulations. 
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Appendix A: Physical quantities for holes 

This Appendix is about the analysis of the Bethe equations in the continuum limit. Here we prove Eq. (|3.20p . which 
gives the variation of a physical quantity A in the presence of a hole A/i in the distribution {A^^''}. It is useful first to 
give the Fourier transform of the momentum and the kernels: 



2fc'(cj) 



27r sinh(7r/2 — 7)0; 



(A.l) 



sinh(7ra;/2) 



^(0)(^) = _27rsinhC 



sinh TTtj 
sinh(7rti;/2) 




27r sinh 27^ 



sinh TTW 



(A.2) 



1 + j(+)(w) 



cosh(7rti;/2) 



(A.3) 



2 sinh(7r/2 — 7)ti; COSI17U; 



2 cosli(7r/2 — 7)(jj cos1i7(jj 
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The hole Xh affects the Lieb equations p.l6p : 

2fc'(A) = [(2^-A-(o))*pW](A)-(ifW*pW)(A) + iA'(")(A-A,), 

2fc'(A) = -(A-(i)*p("))(A) + [(2^-X("))*p(i)](A) + iA^(i)(A-A,,), 
Combming with the ground-state equation, we get: 

p(")(A) + (A) - 2p^(A) = J<+HA - A,) . 

The variation of A is given by: 

dX [p(")(A) - Poo (A)] a(X) + dX [p(i)(A) - poo(A)] a(A) - -a(A^) 

-OO J — OO -^^ 

/OO 1 
dA(pW + - 2poo)(A) a(A) - -a(A;,) 

1 r°° 1 

= -^y_^rfA j(+)(A-A;,)a(A)--a(A,). 
Since J^+^ is even, we get the result p.20p . 



(A.4) 



(A.5) 



(A.6) 



Appendix B: Finite-size corrections 

In this Appendix, we introduce a variant of the Wiener- Hopf method [43], to calculate finite-size corrections to the 
energies from the analysis of the Bethe equations. 



We consider combined magnetic excitations (m'-'^^, m'^^). Since the Bethe integer distributions {if''^^} arc sym- 
metric around zero, the bounds of the integrals in Eq. ((XTB)) arc such that C^"' = C^^^ = B'-^\ Wc can write: 



2fc'(A) = 2^p("' (A) - / dn (/i) A'^"^ (X^fi)- [ dn ifi)K^~^^ (A - ^) 
2A:'(A) 2^p(i'(A)- / dfi p'^°^ {fi)K^^\X - fi) - / p^^) (A - ^i) 

/OO />oo 
Poo(m)^'"'(A - m) - / Poo(m)A"'~'HA - ■ 
-OO ^ — OO 

Applying the convolution by {27t)^^{S + J^^^) to symmetric and antisymmetric combinations, we get: 

(p(o)+pW-2poo)(A) = - / d^ip'^''\^i)J^+\x-^i)- I d^i p^'\^l)j^+\x- ^,) 
(pW-pW)(A) = - / d^Jip('\^i)J^-\x-^l)+ [ df^ p^'\ti)j^-\x- n) . 

J\fi\>B>.o) J\ti\>BW 

Combining again the equations, we get: 

(p(°)-Pco)(A) = - / d/.p(")(M)J(")(A-M)- / d^,p(^H^^)J^-'H^-^^) 

(p«-Pco)(A) = - / dpp^"\f,)J^'\X-^)- f dMpW(^)j(")(A-M). 

J|^|>B(0) J|ai|>S<i) 



(B.l) 



(B.2) 



(B.3) 
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where jW, are defined in (|3.26p . Let us define the symmetric/antisymmetric physical quantities: 

A(±)= / dX p^°\X)a{X) ± dA /i)(A)a(A). (B.4) 



(B.5) 



The variation of A^"^^ with respect to the ground-state value {A^'^'>)ac, can be expressed as: 



d^i (/^) [(J + J(±) ) * a] (m) ± / p(^) (/^) [(5 + ) * a] (/i) 

\fj.\>Bm J\tJ.\>BW 



where we use the fact that J'^^-' and ^ are even. Setting a(A) = —1 or a(A) = e(A), we get the charges and the 
energy: 



m 
N 

fh 



[1 + J(+) (0)] ( / dfi (m) + / rf/i p(^) (Ai) ) , 

\J\H\>BW J|m|>B'" / 



(B.6) 



[1 + J(-) (0)] / (m) - / df, (/i) , (B.7) 



d/^ p(")(Ai)ed(M) + / d^l p(i)(M)e<j(/^) • (B.8) 



To solve the Licb equations (jB.3p . we define the shifted densities: ^'"''(A) ~ + A) for a — 0,1. Neglecting 

the terms from /i < (see [13]), we get the coupled Wiener- Hopf equations: 



9^"\X)+ dMgW(A.)jW(A-A*)+ / dtig^'\fi)j'^-'\X-fi + b) = p^{B(°^ + X) 



5W(A)+/ d^jig^^\p)J^'\X-p-b)+ d^jig^'\p.)j'^°\X-^i) = p^iB^^'^ + X) , 
Jo Jo 



(B.9) 



where b ~ B^^^ — B^^\ After Fourier transform: 



[1 + J(0)(u;)] 5f (c.)+c--^j(-i)(c.) .g^^)(^)+/)(^) = e--s'°'poo(c^) 
e'"'' J(i) (lo) gf {lo) + [l + j'"' g^l^ (uj) + .gi'^ (uj) = c"'"^'" (c^) 



(B.IO) 



We use the factorisations: 



where: 



/ 27r2 r(iw/2) 

G+(tj) = W — ,^ , ' ,^ , , X. 1 , G_(w) = G+(-w) , (B.12) 

+ ^ ^ Y 7r-27 r[(l/2-7/7r)icj] r[l/2 + (7/7r)icj] ' ^' +^ " ^ ' 

H^{u.)= H^{^) = H^{-u) , (B.13) 
r[l/2 -I- (1/2 - 7/7r)icjJ r[l/2 + (7/7r)icjJ 
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and we factorize the 2x2 matrix: 



The matrices G± read: 



1+J 



G+ = ^ 



gi^fc j(l) 1 _^ j(0) 



= G 



1 r^-l 



G± + H± e-'^^Gi - H± 



2 I c'^^(G± -if±) G±+H± 



2 V c'"''(G±' - H^') 



-1 , EJ-l 



G±^ + if. 



Wc can write the system ()B.10|) as: 



Wc multiply by G_: 



GTl I ^+ 



9 



(1) 



V 3 



(1) 



C Poo 



1 



The solution is given in terms of the pole ljq = — i7r/(27) for poo the residue tq = Res(poo7^o)- 



G, 



C Poo [ ^iujb 



--^G,(.)G_(.o) ( 



roCWG_(c^o) 

LUq — LU 



G+(^) 



1 



where C*^"-* = e "^o^' ' . So the magnetic charges are given by: 

m ^ _2roG_(wo).^(o) 
iV woG_(0) 
TO ^ 2roG_(^o) 



The total energy is: 



N 



= 2i^sin27 ^oG-(a;o)G+(-^o) ^^^(o))2 ^ 



27rt; 



27rw 



G_(0)TOy /iJ_(0)m^^ 



TV 



TV 

2 



(B.14) 

(B.15) 
(B.16) 

(B.17) 
(B.18) 



(B.19) 

(B.20) 
(B.21) 



(B.22) 
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Using expressions (jA.3|) for J^^\ we get the critical exponents given in (|3.30p . A similar calculation with C*^"-* 7^ S^"-* 
would give the electric critical exponents. 



Appendix C: Bosonic partition functions and Jacobi's theta functions 
Free boson on a torus 

Let us recall some known results on the free boson theory on a torus [4j. We denote by r the modular ratio of the 
torus, and we write q = e^"^'^. The free boson is defined by the action A and the partition function Zq: 

AM = £ fd'^ IV^P , (C.l) 



Zoig) = exp(-^M) = y^^, (C.2) 

where g is the couphng constant, and 7]{t) is the Dcdekind function: 

00 

,7(T)=gi/24-Q(i_^n)_ (C.3) 
When defects Sip, 6' if are introduced on the boundaries, this defines the partition function Z„i.m', with m,m! integers: 
Zn.M9) = / ,^^2.™ [^'^1 exp(-^M) = Zo{9) exp (^-^^Z^-_!Z^^ . (C.4) 

A Poisson summation of (|C.4p yields: 

Jacobi theta functions 

The Jacobi theta functions are defined as: 

0i(t) = -iE(-l)"(z("+l/2)V2^0, 02(r) = ^q("+i/2)V2^ 



(C.6) 



/2 



nGZ 



They obey the algebraic relations: 



(r) 03(r)04(T) = 2r;(r)^ (C.7) 
V^3(r)04(T) = ei{2T). (C.8) 
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We denote the Jacobi partition functions by: 



?7(r) 



= 2,3,4. 



(C.9) 



The Ising partition functions Z{r, r') are related to the Jacobi ones by: 

Z2 = Z(1,0) + Z(1,1), 

Z3 = Z(0,1) + Z(1,0), 

Z4 - Z(0,1) + Z(1,1), 

Z(0,0) = Z(0,1) + Z(1,0) + Z(1,1), 



(C.IO) 



Appendix D: Partition functions for the staggered models 
Twisted vertex model and Potts model 

Starting from the untwisted partition function Z[g), we can proceed hke in [44j . to construct the twisted partition 
function and the Potts partition function. The partition function where non-contractible loops have a weight n = 2 cos </> 
is given by: 

f(g,0) = 2 L4 ^ +5 +^ +^ H p„^„,'(5) cos(20mAm'), (D.l) 

\ m.in' even ni oven, rn' odd m odd,m' even m,m' odd / 

where m A m' denotes the greatest common factor between m and m' . In particular, for ~ 7r/2, 7r/4, we have: 



Z(ff,V2) - 2|a ^ -S E E 

m,m' even m even,m' odd m odd,/?!' even ni.m' odd y 



Z(g,^/4) - A 



E ^'"'"''(16^)^^ E 



The Q-state Potts partition function has an extra term due to clusters with cross geometry [44] : 



where 



ZpottsiQ) = Z{g, neo) + ^{Q - l)Z{g, n/2) , 



/— - TT TT — 27 

VQ = 2cos7, 0<7<-, g=^— 



7 

eo = - . 

TT 



(D.2) 
(D.3) 

(D.4) 
(D.5) 



Particular values of Q 
• The case Q = 2. 

This provides a good check of the result (|4.14p . since the Potts model arising from the staggered vertex model 
is equivalent, on the lattice, to the usual critical Ising model. Using (jD.4[) : 



^Potts(Q — 2) 



(A-B) J2 +(^-^) E +iA-D) Y 



m even.m' odd 



m odd.m' even 



m.jn' odd 



^m,m'(V4)- (D.6) 
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Now the sums on m, m' can be expressed in terms of the Zi^: 

m even,m' odd ra odd,m' even m.m' odd 

(D.7) 

We obtained the first identity by using (jC.7p . and the two others by expanding the square of the left-hand sides. 
Combining pjp with (gjl]) and (|C.8|) . we get: 

^Potts(Q = 2) = 2^^^ + Z3 + Z4) Rising , (D.8) 

so we correctly find the Ising partition function. 
• The case (5=1. 

This case is a priori a bit intriguing. The partition function of the Potts model is then a trivial object (since 
there is only one state available for the whole lattice), while the general formulas for the central charge give 
in this particular case c = —2 (Q = 1 so 7 = 7r/3, g = 1/6, eo = 1/3). This discrepancy occurs for the same 
reason as in the Berker-Kadanoff phase [5] : the level of the transfer matrix corresponding to a trivial partition 
function (and hence, formally, c = 0) is very high in the spectrum, while the level gcnerically dominating the 
thermodynamics (but which disappears right at Q = 1 hy quantum group truncation) corresponds to c = — 2 
(this means the free energy is a discontinuous function of Q or of the boundary conditions [3]). 

Let us now see the mechanism in more details. The ground-state energy of our system in the untwisted case 
is twice the ground state energy of the antifcrromagnetic XXZ model with Aq = — cos 2j. In the case Q — 1 
we have Aq = 5. The antifcrromagnetic XXZ model with this value of the anisotropy is related with the Potts 
model at Q = 1 on the 'non-physical self-dual line' [32]. Recall that, meanwhile, the Potts model on the usual 
self-dual line is related to the antifcrromagnetic XXZ chain at 2 A = —\fQ. so A = — ^ in the case Q = \. 

Now wc know that the energies of the antifcrromagnetic XXZ at A = — i are minus the energies of the anti- 
ferromagnetic XXZ at A = i (this is the general mapping between and — i?_A). The ground-state energy 
of the antifcrromagnetic XXZ at A = ^ is the same, per unit length in the thermodynamic limit, as the one of 
the twisted antifcrromagnetic XXZ, i.e. the ground-state energy of the percolation problem, i.e. Eq = Q in the 
proper normalisation. We thus conclude that the eigenvalue 'corresponding to Z = 1' in our spectrum is the 
most excited among the .subset of symmetric states. 

It is useful to see this mechanism at the level of partition functions as well. Start from (jD.l[) and set Q = l,4> = 
7r/3. Then there is no contribution from the cross-geometry clusters. Since cos27r/3 = cos47r/3 = — i, we have 

Zm,m'{g) COS(20 mAm') = - ^ ^m,m' ^ 2 ^ ^m,m' • 

m.m' m—3p,m'—3p' m.m' 

Moreover, m — 3p is odd (resp. even) iff p is odd (resp. even). Finally, ^3p,3p'(g) = Zpy{9g)/3. So we can 
rewrite 

^Potts(Q = 1) = U J2 +^ +^ J2 +^ J2 [^™,m'(9.9)-Z„^,„'(5)] . (D.9) 

\ 7n,m' even m even,m' odd m odd,m' even m,7n' odd / 
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We can recombine terms using expressions for the A, B, C, D in terms of the Z^. We find 

^Potts(Q = l) = ^{A~B^C^D) + iB-D) [Zm.m'{3/2)-Zm,ra'{l/6)] 

m cvcn,m' 

+ (C-D) [^m,m'(3/2)-Z™,™'(l/6)] , (D.IO) 

m,m' even 

where we have speciahzed to g = 1/6 and used Eulcr's identity: 

Y ^,n,m' (3/2) (1/6) = i[Z,(6)-Z,((2/3)] = l. (D.ll) 

ni^m' even 

Now we have 

Y [^™,rn'(3/2)-Z,„,,„,(l/6)] = [^em(3/2)-Ze,„(l/6)] (D.12) 

771 even,m' 7/1 cvcn,e 

and 

Y [^m,rn'(3/2)-Z„,,„Kl/6)] = ^ Y [^e,n (3/2) - (1/6)] . (D.13) 

ra.m' even m,e half — integer 

Both terms can be shown to vanish exactly. We conclude that 

^Potts(Q = 1) = 0. (D.14) 

This means there are exact cancellations among states in the low-energy spectrum: the (unique) state that 
would correspond to the trivial partition function is very highly excited and does not contribute to the conformal 
partition function (at c = — 2 in this case). 
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